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I. INTRODUCTION 

In the last decade, the increasing energies of the collision experiments have leaded to 
the detection of various new and "exotic" mesons as, for example, the K(4140) (reported by 
CDF Collaboration as a narrow near-threshold structure in the J/ip<p mass spectrum PQ) and 
the X(4350) (reported by Belle Collaboration also as another narrow structure in the J/ip(j) 
mass spectrum at 4.35 GeV [2]). By exotic, we mean that it is believed that these mesons are 
beyond the usual quark-model description as qq pairs. Since most of their quantum numbers 
have not been measured yet, there are many uncertainties about their constitution. Various 
models started to appear attempting to explain the masses and the observed decay modes 
of these mesons. Some of these models consider these exotic mesons as quark-gluon hybrids 
(qqg), tetraquark states (qqqq), molecular states of two ordinary mesons, glueballs, states 
with exotic quantum numbers and many others [3H7]. Some of these hypothesis have been 
studied in many works using the QCD Sum Rules (QCDSR) approach regarding masses, 
branching ratios and decay constants, obtaining different degrees of success depending on 
the choice of the interpolating currents [BHT2]. 

Following the development of previous works of our group ( [T3HT5] and references therein, 
just to name a few), in this work we use the QCDSR formalism to obtain the coupling 
constants of the meson vertices J/ipD s D s and <pD s D s . These vertices may appear in the 
decay processes y(4140) — > J/ijjcj) and X(4350) — > J/ip(f>, depending on the model used to 
represent the y(4140) and X(4350) mesons. We call the attention to the presence of these 
two vertices in the works of Ref. [5], where the observed decay F(4140) — > J/ip(p is studied 
with the intermediate decay y(4140) — > DgDi*^ — > J/ip<p, and Ref. [7], where the observed 
decay X(4350) — > J/ip<p is analyzed using the intermediate process X(4350) D^dP 
J/i/j(j). It is worth to mention that a more precise knowledge of these coupling constants 
and form factors certainly will improve the understanding of the fundamental constitution 
of these new observed mesons. 

II. FORMALISM 

The starting point for a three meson vertex in QCDSR is the three point correlation 
function [15] . In this work, we are interested in two vertices of the type vector-pseudoscalar- 
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pseudoscalar (VPP), namely the J/ipD s D s and <pD s D s vertices. Therefore, we have two 
different correlation functions for each vertex, one with a vector meson off-shell and another 
with the pseudo-scalar meson (D s ) off-shell: 

T^W) = / (0'\T{ 3 ^(x)j^\y)j^\0)}\0')e^e-^d 4 xd 4 y (f) 



rf* W) = / V\T{j?'(x)j?' \y)j^(0)}\0')e^e-^d 4 xd 4 y n = 1, 2 (2) 



where V\ = <f> and V2 = J ftp, being q = p' — p the transferred momentum. 

According to the QCDSR, we can calculate these correlation functions in two different 
ways: with hadron degrees of freedom (called the phenomenological side) or with quark 
degrees of freedom (called the OPE side). Both representations are equivalent and can be 
equated invoking the quark-hadron duality, after the application of a double Borel transform 
on each side. Then, the form factors and coupling constants can be extracted. 

A. The phenomenological side 

In order to accomplish the calculation of the phenomenological side, it is necessary to 
know the effective Lagrangians of the interaction which, for the vertices V n D s D s , are Ej: 

C VnDs D B = {-l) n ig Vn D 3 D s V:{Dtd a D- - d a D}Dj), n = 1, 2 (3) 

From these Lagrangians, we can obtain the vertex of the hadronic process. In the case of 
V n off-shell we have: 

(D s (p)v n (q)\D s (p)} = (-i) n i S fi$. D .(<t t )f{q)(p a +P' a ) (4) 

We also make use of the hadronic matrix elements 



(0\j^\D s (p')) = (D s (p)\j°°\0) = f Ds -^- (5) 
(V n (q)\j^\0) = f Vn m Vn e;(q) (6) 

In the case of D s off-shell the obtained vertex is: 

(D s (q)V n (p)\D s (p')) = (-l)^i D M 2 )e a (p)(2p' a - p a ) (7) 

and the corresponding usual matrix elements are: 
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(0\j?°\D s (p')) = (2?.(g)L?f -|0> = fn.—*- (8) 

/ / Lf~ \~ 1 1 be 



(K(p)b;lO) = fv n mv n e;(p) (9) 
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In the formulas above, fy n and fo s are the decay constants of the vector (n = 1,2) and 
D s mesons respectively, e M ( a ) is the polarization vector, D (g 2 ) and gy°D d (l 2 ) are the 
form factors of the vertices, which correspond to V n or D s messon off-shell respectivelly, m s , 
m c are the strange and charm quark masses and my n is the mass of the vector meson. 

In order to improve the matching between the phenomenological and QCD sides, it is 
convenient to make the change of variables p 2 ->■ -P 2 , p' 2 -> -P' 2 and q 2 ->■ -Q 2 . After 
that, we can obtain the correlation functions for the phenomenological side in the form: 



r phen(V n ) _ ^YVvikDs {Q 2 )fl a fv^ D m Vn ( Pfl + 

*■ II / \ n / T,n Ow^n. ON / t-..o n \ ~T~ (1.7 . 



(m c + m s ) 2 (P 2 + m 2 D ){Q 2 + m 2 v )(P' 2 + m 2 



(10) 



h (D , {-^) n 9tinM)flsfy n rn%my n (^=^ - 2p' 

-rphen(D s ) V Yn Zj_ \U r 11 

" {m c + m s ) 2 {P 2 + m 2 v J{Q 2 + m 2 D J{P' 2 + m 2 Ds ) 1 j 

where h.r. stands for the contributions of higher resonances and continuum states of the 
involved mesons and <?y£> d (Q 2 ) * s the f° rm factor of the V n D s D s vertex with meson M 
(M = V n ,D s ) off-shell. We remind that the limit lim d (Q 2 ) * s our definition for 

the coupling constant gy n D s D s of the vertex [TBI . 



B. The OPE side 



The OPE side is obtained using the interpolating currents written in terms of the quark 
fields in Eqs. (jl| and (2 ). The meson interpolating currents used in this work are j J J^ = cj^c, 

By construction, the OPE side is given by an expansion known as Wilson's Operator 
Product Expansion. In the form factor calculation, this expansion exhibits a rapid conver- 
gence and can be truncated after a few terms. In this work, we consider contributions from 
two kinds of terms in the OPE side: the perturbative (the dominant one) and the quark 
condensates. In previous works of our group |T3l [T7] , it was shown that the gluon condensate 
term has a negligible contribution. For this reason, it has not been included in the present 
work. By using dispersion relations, the correlation function for a given meson M off-shell 
can be written in the following form: 

rOPBM M= l / P* ^'V dsdu + rMM (12) 

^ Jo Jo {s - p 2 ){u - p' 2 ) 
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where the first and second terms are the perturbative and quark-condensate contributions 
respectively. The spectral density p^ rt ^ M ' (s,u,t) is the double discontinuity of the pertur- 
bative term of the OPE, obtained by the use of the Cutkosky's rules and s = p 2 , u = p' 2 , 
and t = q 2 are the Mandelstam variables of the vertex. 

Invoking Lorentz symmetry, the spectral density for a VPP vertex can be parametrized 

as: 



pPert(M) { 



s, u, t) 



where A = (u + s — t) 2 — 4ms, and Fp and F^ vl> are invariant amplitudes. For the cases 
studied in this work, these invariant amplitudes can be written as: 



Fi M \s,u,t)p„ + F( M \s,uM 



(13) 



(M) 



= {-l) n [u{A - 1) - A(t - s) - (2 A - \){m s - 



-l) n [s{B - 1) + B(u - t) + (1 - 2B){m s 



m r 



= u{A + (-1)™) + A(t - s) - ((-l) n + 2A){m s - m c ) 2 



-s(B + (-l) n ) + B(u + t)- 2B{m s 



m. 



(14) 
(15) 
(16) 
(17) 



where 
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k p' Q \k\cos6 
" \p'\Vs 



\k\cOS9 j: 

B = — =; — cos 
\p'\ 



\k\ = y k 2 , — \2 — n — e\m 2 s — |1 — n + e\m 2 k 
s + u — t , ->., \/A 



_ 2p' k -u+ (-lY n+t \m 2 s - m 2 c ) 

2\j}\\k\ 
s + e(-l)( n+6 \m 2 c -m 2 s ) 



\p'\ 



2^~s * 1 2^7s 

and e = 0(1) for D s (V n ) off-shell. 

Diagrammatically, the perturbative contributions are shown in Fig. [T] The quark con- 



densates contributing to Eq. (12) are 



r (s->W = _| n + e _ 2 |(s a ). 



[mcPn + ((e-l)m s + em c )p{J 

1 



(18) 



(p 2 — (1 — e)m 2 s — em 2 c )(p' 2 — m\ 

where (ss) = (0.8 ± 0.2) (qq) = (0.8 ± 0.2) (0.23 ± 0.03) 3 GeV 3 [Sj . The quark condensate 
contributions are represented by the diagrams of Fig. [2j 
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FIG. 1. Perturbative contributions to vertex J/ipD s D s with J/tp off-shell (panel (a)) and D s 
off-shell (panel (b)), and to vertex 4>D S D S with cf) off-shell (panel (c)) and D s off-shell (panel (d)). 




FIG. 2. Quark condensates contributing to vertices J/ipD s D s (panel (a)) and <pD s D s (panel (b)). 



C. The sum rules 



To obtain the sum rules, two more steps are necessary. The first one is to apply a double 



Borel transform [181 [TU] to both sides of the QCDSR in Eqs. (|To|), (jTTJ and Q, which 
involves the following transformation of variables: P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — > M' 2 , 
where M and M' are the Borel masses. After that, we can equate both sides of the QCDSR 
invoking the quark-hadron duality. Because of the double Borel transform, quark condensate 
contributions to vertices V n D s D s are negleted when the lighter meson is off-shell. This has 



been taken into account when writing Eq. (18). 



The second step is to eliminate the h.r. terms appearing in the phenomenological side 



in Eqs. (10) and (11). This is achieved by the introduction of the continuum threshold 
parameters So and uq. These parameters satisfies the following relations: mf < Sq < m! 2 
and ml < Uq < m' 2 , where m; and m are the masses of the incoming and outcoming mesons 
respectively and m' is the mass of the first excited state of these mesons. Taking advantage 
of the quark-hadron duality, the h.r. terms can be ride away from the QCDSR. 

After performing these steps, we can match the phenomenological and QCD sides, ob- 
taining the QCDSR expressions for the form factors. For example, in the case of the p^ 
structure for the V n meson off-shell and the p'^ structure for the D s meson off-shell, the 
obtained form factors are: 

2 2 

-A P P -yfFi Vn) e-^e-^dsdu- In- l)m c (ss)e"^e"^ 

(Vn) r n 2\ _ 8tt 2 Js inf Ju inf p V ' CX ' /-.q\ 

(m c +m s ) 2 (Q 2 +m^) e " e 

2 2 

-A f S ° f U ° J rrF^ s) e~^e~^dsdu + (2 - n)m s (ss)e~^ 

JDs) (n 2\ 8tt 2 Js inf Ju inf yg P V J s\ / 

9VnD.D.W ) - (_ 1} n 2/ 2 ^ ^ ^ ^ ^ /M2 /M , 2 {^) 

(m c +m s ) 2 (Q 2 +m 2 Dg ) 

remembering that n — 1 for V± — <p and n = 2 for Vi = J/ip. 



III. RESULTS AND DISCUSSION 



Eqs. (19) and (20) show the two different form factors for each vertex. In order to 
minimize the uncertainties when extrapolating the QCDSR results, it is imposed that both 
form factors lead to the same coupling constant [16]. This condition was used to reduce 
the errors of finding the Borel masses and the continuum thresholds, which should gives a 
good plateau in order to obtain a physical results. Form factors should not depend on these 
quantities. A good plateau is the result of a sum rule which gives a good stability. The 
plateau is often called "window of stability" . 

Besides the values of the Borel masses and the continuum thresholds, we need to know 
the values of the hadron masses and decay constants. Their values are given in Table |TJ 

2 

The Borel masses, which respect the relation M' 2 = M 2 , can assume any value within 
the Borel window. This window is obtaining after imposing that the pole contribution must 
be bigger that the continuum contribution by 50% and the quark condensate term contribute 
only with 30 % of the perturbative term. We have worked with the mean values of the form 
factors within a Borel window. The procedure is briefly sketched as follows: the mean value 



TABLE I. Parameters used in this work, taken from Ref. [20J. 



s c (p D s J ftp 



m (GeV) 


0.101 


1.27 


1.019 


1.96847 


3.096916 


/ (MeV) 






229 


257 


416 



of the form factor is calculated in the Borel window for each value of Q 2 used. The standard 
deviation is used to automatize the analysis of the stability of the form factor with respect 
to the Borel masses and continuum threshold parameters. In this way, it is guarantee a good 
stability in the Borel window and in the whole Q 2 interval. 

sq and «o are defining as Sq = (rrii + Aj) 2 and Uq = (m + A ) 2 , where the quantities 
Aj and A„ have been obtained imposing the most stable sum rule. In order to include the 
pole and to exclude the h.r. contributions for the cases of 0, D s , and J /if) mesons off-shell, 
the values for Aj/^, A^,, and Ad s cannot be far from the experimental value of the distance 
between the pole and the first excited state (201 EI] • Our analysis has found that the best 
values are A^ = 0.6 GeV, A Ds = 0.6 GeV and Aj/^ = 0.5 GeV, which leads to a remarkably 
stable Borel windows, as can be seen in Fig. [3] for the J/ipD s D s case. 

About the choice of the Dirac structures used in the calculations, it is worthy to say that 
in principle, and if the whole OPE series could be summed up, both structures, p^ and p'^ in 



Eq. (13) would lead to a valid sum rule. In a real calculation however, the OPE series has to 
be truncated at some order, and approximations are necessary to deal with the h.r. terms. 
Thus it is not always possible to use both structures any more. Regarding the structure p M 
for D s off-shell, it does not lead to a coupling constant compatible with the <f> off-shell case 
in the vertex <pD s D s . For the others three form factors, namely <7jmp s £> s (Q 2 ), 9md s d s (Q 2 ) 
and gsD s D 3 {Q 2 )i ^ ne structure p^ leads to results very close to the obtained using the p'^ 
structure. This was taken into account when calculating the uncertainties in this work. 



A. Results for the J/tpD s D s and <fiD s D s vertices 



In Table 



II 



we present the Q 2 and Borel windows found for each form factor 9v^d s d s (Q 2 )^ 
together with its parametrization and the respective coupling constant gv n D s D a with its 
associated error a, calculated using the method explained in subsection B. 
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M 2 (GeV 2 ) 



FIG. 3. Form factor stability for the J/^D S D S vertex at Q 2 = 1 GeV 2 



TABLE II. Parametrization of the form factors and numerical results for the coupling constant of 
this work. The calculation of a is explained in subsection B. 



Quantity 



Vertex 



J/ipD s D s 
J /ip off-shell D s off-shell 



(f>D s D s 

off-shell D s off-shell 



Q 2 (GeV 2 ) 
M 2 (GeV 2 ) 
a {M) (O 2 ) 
A 

B (GeV 2 ) 

(M) i 
9v n D s D s ± ° 



[0.5, 2.0] 
[3.1, 4.6] 

3.21 
15.52 



[1.0, 5.0] 
[3.0, 4.7] 

2.503 
4.425 



5.96 



-0.97 
-0.91 



6.01 



+0.52 
-0.43 



[1.0, 2.5] 

[1.2, 2.0] 

Ae -Q 2 /B 

1.091 
1.889 



[0.1, 1.6] 
[0.6, 1.3] 



B+Q 7 
16.25 GeV 2 

13.3 



1.89 



+0.18 
-0.14 



1.73 



+0.09 
-0.09 



The Borel windows presented in Table |TT] (M 2 row) satisfies the already mentioned con- 
ditions regarding the pole and continuum contributions, as can be seen in Fig. [4] for the 
J/ipD s D s vertex. For the cj)D s D s case we obtained similar results. These Borel windows 
also satisfies the dominance relations between the perturbative over the condensate contri- 
butions by at least 70% of the total contribution for both vertices. 

The form factors obtained for the J ji\)D s D s vertex, with J/ijj and D s off-shell, were both 
well adjusted by exponential curves (Fig. pj^a)). For the vertex cj)D s D s , the form factor for 
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FIG. 4. Pole and continuum contributions for the J/tp off-shell (panel (a)) and for D s off-shell 
(panel (b)), both for the J/4>D S D S vertex at Q 2 = 1 GeV 2 . 

the case D s off-shell was fitted by a monopolar curve while for off-shell it was adjusted by 
an exponential curve (Fig. [5](b)). 




-10 -8 -6 -4 -2 2 4 6 -4 -3 -2 -1 1 2 3 

Q 2 (GeV 2 ) Q 2 (GeV 2 ) 
(a) (b) 



FIG. 5. Form factors for the J/ipD s D s vertex (panel (a)) and for the (f)D s D s vertex (panel (b)). 



B. Error analysis 

The vertex J/i/jD s D s has similarities with the vertex cf)D s D s . However, they are physically 
different in such a way that it is expected different values for their coupling constants. 
Nevertheless, the two coupling constants of the same vertex should have the same value, 
independently of which meson is off-shell. The coupling constants presented in Table [IT] 
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have different values among them even when we look for the two coupling constants of 
the same vertex. But when the errors are taken into account, we can see that they are 
compatible, as suggested by Fig. [5j where the error bars stands for the uncertainties of the 
coupling constants in Table [TT} 

In works [HI [22], we estimated the errors by studying the behavior of the coupling 
constant when each one of the parameters involved in the calculation were varied individu- 
ally between their upper and lower limits. Experimental parameters have their own errors 



(showed in Table III), while for the QCDSR parameters, we varied the thresholds A s and 
A u in ±0.1 GeV, the momentum Q 2 in ±20% and for error due to the Borel mass M 2 we 
used the standard deviation of the average value of the form factor within the Borel window. 
Then, we calculated the mean value of all these contributions and their standard deviations. 
The same procedure was adopted in this work in other to estimate the error bars. In Ta- 



ble III , we can see how the variation of each parameter affects the final value of the coupling 
constants. 



TABLE III. Percentage deviation of the coupling constants related with each parameter. 











Deviation 


(%) 




Parameter 














f Da = 257.5 ±6.1(MeV) 


3.88 


3.87 


3.88 


3.86 


f m = 416 ± 6(MeV) 


1.18 


1.18 






ft = 229 ± 3(MeV) 








1.07 


1.07 


m c = 1.271^9 ( GeV ) 


14.7 


18.9 


15.4 


11.1 


m s = lOljJ? (MeV) 


1.25 


3.10 


5.43 


3.37 


M 2 (GeV 2 jj 






6.46 


7.70 


6.56 


0.72 


A s ±0.1(GeV) 






32.5 


2.70 


9.56 


3.09 


A„±0.1(GeV) 






26.7 


2.08 


12.0 


0.46 


Q 2 ± 20% (GeV 2 ) 


a 




3.38 


2.34 


2.91 


1.54 


(ss) = (0.8 ± 0.2)(0.23 ± 0.03) 3 (GeV 3 ) 


7.70 






5.00 


Fitting parameters (A and B) 


0.94 


0.11 


0.15 


3.12 



The intervals for these quantities are those of Table UJ\ 
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It follows from Table III that the variation of most of the parameters (namely, the decay 
constants, the Borel and the strange quark masses for example) have little impact on the 
value of the coupling constants. We have found however that the error in the charm mass is 
the one with the biggest propagation over the final value of the coupling constants. The fact 
that the sensitivity related to the Borel mass is small was already expected, as we obtained 
a very stable Borel windows for both vertices (Fig. [3]). 

All the results showed until this point are for the structure p^ in the case V n off-shell 
and for the structure p'^ in the case D s off-shell. As said earlier, we can also work with 
other (good) structures in order to obtain the final value of gj/^D s D a and g^D a D s and their 
uncertainties. In the case of J/ipD s D s , we used the coupling constant obtained from the 
study of the structure p^ in the D s off-shell case, which value is: 

9f,i DsDs = 6.561^? (21) 
and also the structure p'^ for the J/ip off-shell case, which value is: 

9 { j%%d s = 6.10j££ (22) 

For the off-shell case in the vertex <j)D s D s we used also the result obtained from the p'^ 
structure, which reads: 

,(*) - i on+o.i7 



$U = L90 -oi^ (23) 



IV. CONCLUSION 



Using the results presented in Table |ll| and in Eqs. (21), (22) and (23), we taken the 
mean value of the calculated coupling constants, obtaining the following final results for the 
coupling constants of the two vertices studied in this work: 

ffW. = 6.20±?;?5 (24) 
9<j>D a D s = 1.85^23 (25) 

We call the attention again to the fact that these results were obtained from a "good" 
QCDSR, which means a good pole/continuum dominance, a perturbative contribution 
greater than the condensate contribution and a very stable Borel window for the whole Q 2 
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interval studied. Also, the coupling constants with the D s meson off-shell were compatible 
with the coupling constants obtained for the vector meson off-shell case in both vertices. 
Regarding the form factors, we observe the same behavior found in our previous works: the 
form factor is harder when the heavier meson is off-shell. The errors are of the order of 
10%-20%, what are the expected values coming from previous QCDSR works. 



The results of Eqs. (24) and (25) can be compared with the ones coming from other 
calculations. For example, the result for the g^DsDs of this work can be directly compared 
with the result coming from the Light Cone QCDSR (LCQCDSR) formalism, as it is shown 



in Table IV, concluding that they are in agreement. We can use also previous QCDSR 



TABLE IV. Meson coupling constants obtained from other QCDSR calculations. 



QCDSR [IT] 


QCDSR [L3] 


LCQCDSR [23] 


dj/ipDD = 5.8 ± 0.8 


gjipD-D* = 6.2 ± 0.9 


94>d s d 3 = L45 ±0.34 



results for other vertices to compare with our gj/^D s D s value. For example, invoking the 
(broken) SU(3) symmetry [24J, it is expected the relation gj/^D s D s = Qj/^dd- Using the 



value for gj/^DD from Table IV, we can see that this relation is sustained. Still inside the 
SU(3) scheme, the relation gj/^D a D a = gj/*(>D*D* should be valid [151 121] • Comparing our 



result for g,j/^D a D s with the one for gj/^,D*D* from Table IV, we obtain again compatible 
results. The gj/^D a D s = gj/^DD SU(3) relation is violated by the order of 8%, the same 
value found between gj 1 po*D* and gj/^DD [15] . which means that, in spite of the huge mass 
difference between the involved mesons, SU(3) is a reasonable symmetry for this vertices 
when using QCDSR. 
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